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Employing first principles density functional calculations, together with solution of the low-energy, 
model Hamiltonian constructed in a first principles manner, we explored the origin of magnetism 
and Tc trend in Cr-based double perovskite series, Sr2CrB'06 (B'=W/Re/Os). Our study shows 
that the apparently puzzling Tc trend in Sr2CrB'06 (B'=W/Re/Os) series can be understood in 
terms of the interplay of the hybridization driven mechanism and the super-exchange mechanism. 

PACS numbers: 71.20.-b, 71.20.Be, 75.50.-y 



I. INTRODUCTION 

Oxides with high magnetic transition temperature (T^) 
are important for technological advancement. A much 
discussed family of compounds in this connection are 
double perovskites, with a general formula A2BB'06, 
A=alkaline/rare-earth metals and B/B'=transition met- 
als. The observatioEpI of large magnetoresistance with 
a fairly high ferromegnetic (FM) Tc of about 410 K in 
Sr2FeMo06 (SFMO) and its unusual origin of magnetism 
brought this family of compounds into the forefront of 
activity. Since then the question has been, can Tc be 
boosted even further and if so, is there a systematic 
trend that can be observed and understood. Attempts 
have been made to dope SFMO with La to increase the 
Tc- Though some partial success has been achievecP fol- 
lowing this path, Tc was found to be boosted much more 
efficiently by moving to different choices of B and B' ions. 
The microscopic understanding of this increase, however, 
has not been achieved. A collection of measured transi- 
tion temperatures seems to bear some correlation with 
the number of valence electrons among different double 
perovskite compounds (see Fig.l in RefO]). Focusing onto 
Cr-based series, namely Sr2CrB'06 (B' = W/Re/Os), the 
family with spectacularly high Tc, the measured Tc shows 
a rapid increase as one moves from SraCrWOe (SCWOj^l 
with Tc ~ 450 K to SraCrReOe (SCRO)^ with Tc w 620 
K to SrzCrOsOe (SCOO)^- with Tc w 725 K. This rapid 
increase is curio us considering the hybridization driven 
(HD) mechanisnPElof magnetism accepted so far for dou- 
ble perovskites, proposed in the context of SFMO, which 
would predict decrease rather than the increase.™ Fur- 
thermore, the compound with the highest Tc in this se- 
ries, SCOO, is found to be an insulator making the issue 
even more interesting. Also, a distinction between oth- 
ers like SFMO and La doped SFMO and this series is 
a possible role of spin-orbit coupling (SOC) due to the 
presence of 5d elements. 

We have used density functional theory (DFT) based 
calculations together with exact diagonalization of Cr- 
B' model Hamiltonian constructed in a first-principles- 
derived Wannier function basis, in order to probe this 



issue. Though few DFT based studiei^ exist for indi- 
viduals of these compounds, to best of our knowledge, 
theoretically no comprehensive study exists to address 
the origin of growing trend in Tc within the series. 



II. METHODOLOGY 

The DFT calculations were carried out using the 
plane wave pseudopotential method implemented within 
Vienna Ab-initio Simulation Package(VASP)."'^^ The 
exchange-correlation functionals were approximated by 
generalized gradient approximation (GGA).'i2l We have 
used projected augmented wave (PAW) potentialPI and 
the kinetic energy cut-off for expansion of wavefunctions 
used was 450 eV. Reciprocal space integrations have been 
carried out with a k-space mesh of 6x6x6. For ex- 
traction of a few-band, tight-binding Hamiltonian out of 
full GGA calculation, we have carried out muffin-tin or- 
bital (MTO) based NMTO-downfolding calculations .231 
The reliability of the calculations in two basis sets has 
been cross-checked. 



III. RESULTS 

A. Examination of Basic Electronic Structure 

In order to unravel the origin of magnetism in 
Sr2CrB'06 series, let us first critically examine the elec- 
tronic density of states (DOS) of these compounds. Left 
panel of Fig.l shows the DOS, as obtained in spin- 
polarized DFT calculations within GGA. The states close 
to Fermi level {Ep) are dominated by Cr and B' d states 
hybridized with O p states, while the O p dominated 
states separated from Cr and B' d dominated states oc- 
cupy the energy range far below Ep and Sr s and d dom- 
inated state remain far above Ep. The d states of Cr 
and the B' ions are exchange split as well as crystal field 
split. The empty B' t2g states in the up-spin channel ap- 
pear in between the crystal field split Cr t2g and eg states, 
gaped from Ep while the B' states in the down spin chan- 
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Figure 1: (Color online) Left Panels: GGA DOS projected 
onto Cr d (black solid lines), B' d (Red/Grey solid lines) and 
O p (shaded area). Zero of the energy is set at Ep. The num- 
bers within the boxes indicate the Cr t2g contribution in the 
bands crossing Ep, with respect to that of B'. Right Panel: 
Calculated net magnetic moment and magnetic moments at 
Cr, B' and O sites. The numbers within the bracket denote 
the result of GGA+SO calculations, the first entry being the 
spin moment and the second entry being the orbital moment. 
From top to bottom, the plots correspond to SCWO, SCRO 
and SCOO respectively. 
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Figure 2: (Color online) The energy level diagram (upper 
panel) and massively downfolded Wannier functions (lower 
panels) for Sr2CrB'06 series. For Wannier function plots, 
constant value surfaces have been plotted with two oppositely 
signed lobes colored differently. From left to right in lower 
panel, the plots correspond to SCWO, SCRO and SCOO re- 
spectively. The numbers in the energy level diagram are in 
unit of eV. 



B ' site following the dehybridization effect between Cr 
and B'. 



nel hybridized with Cr t2g states, either cross the Fermi 
level as in the case of W and Re compounds, or remain 
completely occupied, as in case of the Os compound. It is 
rather intriguing to notice that the hybridization between 
Cr t2g and B' t2g in the down spin channel, progressively 
gets weakened in moving from W to Re to Os compound. 
This may be appreciated by considering the Cr contribu- 
tion measured with respect to B' contribution for the 
states in the down spin channel close to Ep. For W, Re 
and Os compounds, it is found to be 66%, 16 % and 5% 
respectively. This is caused by the gradual moving down 
of the B ' energy level, in moving from left to right of 
the periodic table across the same row (W Re ^ Os), 
reflecting an increase in the ionic potential experienced 
by the 5d electrons with an increase in nuclear charge. 
As discussed later, the hopping interaction connecting Cr 
and B' t2g states on the other hand remains similar across 
the series. Analysing the DFT calculated magnetic mo- 
ments, presented in the right panel of Fig.l, we find the 
second interesting observation that the magnetic moment 
per electron at B' site, defined as m/d, where m is the 
calculated moment and d is the valence count at B ',^51 
keeps growing from W to Re to Os. The increase in m/d 
becomes even more evident taking into account the mag- 
netic moment at O site, which is small and point to Cr 
moment for SCWO, small (large) and point to Re (Os) 
moment for SCRO (SCOO). This prompts us to conclude 
that there is a growing intrinsic moment that develops at 



B. NMTO-downfoIding calculations 

In order to further analyse the findings of the elec- 
tronic structure calculations, we have carried out NMTO 
downfolding calculations, which is engineered to define 
energy-selected, effective Wannier functions by integrat- 
ing out degrees of freedom that are not-of-interest (down- 
folding). As a first step we downfolded O p, Sr as well 
as Cr and B' Cg degrees of freedom. This defines an ef- 
fective basis consisting of Cr t2g and B' t2g states. In 
the second step, we applied massive downfolding, keep- 
ing only B' t2g degrees of freedom active and downfold- 
ing all the rest including Cr t2g degrees of freedom. On 
site matrix elements of the real space Hamiltonian de- 
fined in the Cr t2g-^' t2g basis and the massively down- 
folded basis give the energy level positions before and 
after switching on the hybridization between Cr and B 
' states, respectively. Fig. 2 summarizes the results for 
W, Re and Os compounds. The energy levels in the left 
and right half of the upper panel, demarked with verti- 
cal grey line, depict the energy level positions obtained 
out of downfolded Cr t2g-B' t2g basis and the massively 
downfolded basis, respectively. The lower panels exhibit 
the plots of the one of the t2g (xy) Wannier functions cor- 
responding to massively downfolded Hamiltonian in the 
down spin channel. Examination of Fig. 2 brings out two 
aspects: Firstly, the progressive dehybridization effect, 
as discussed in the context of DOS plots, is evident in 
the plots of Wannier functions. The central parts of the 



3 



Wannier functions are shaped according to B' xy sym- 
metry and the tails of the Wannier functions sitting at 
neighboring sites are shaped according to O p and Cr 
t2g symmetry. The tails reflecting the hybridization be- 
tween the Cr t2g and B' t2g-0 weaken as one moves from 
W to Re to Os compound. As a consequence, the ratio of 
renormalized spin splitting to that of the bare splitting 
at B' site reduces drastically from W (14.5) to Re (3.3) to 
Os (2.2). Secondly, considering the level splitting at B' 
site before switching on the hybridization, we find while 
the splitting at W is negligibly small confirming the non- 
magnetic character of the B' site, those for Re and Os 
are found to be « 0.31 eV and « 0.53 eV, respectively. 
These values are significantly larger compared to what 
one would have expected considering (P valence in Re 
and valence in Os compared to valence in case of 
W with 0.06 eV splitting, which would have given rise to 
splittings of 0.12 eV and 0.18 eV respectively. This con- 
firms the presence of a growing intrinsic, local moment 
at B' site as one moves from W to Re to Os, driven by 
the dehybridization effect. The magnetism in Cr-B' se- 
ries, therefore, needs to be understood as an interplay of 
two mechanisms: HD mechanism as operative in SFMO 
which causes renormalized, negative spin splitting within 
B' states that appear in between the exchange split Cr 
states, and the superexchange (SE) between the mo- 
ment at Cr site and the intrinsic moment at B' site, which 
would align the moments at Cr and B' sites antiparallely. 
For W, the intrinsic moment being negligible, the mag- 
netism is entirely driven by HD mechanism, while for the 
other extreme of Os, the hybridization effect is weak, SE 
having a rather large contribution. The presence of such 
intrinsic moment at 5d site is counterintuitive at a first 
glance. Comparing the situation, with the double per- 
ovskite Sr2ScRe06, for which the magnetism has been 
recently investigated,'^^ Re was found to possess a rather 
small intrinsic moment of size 0.013 iib- Our electronic 
structure calculations carried out for Sr2ScRe06 find also 
a similarly small moment (0.03 /is). The unusual local- 
ized aspect of Re or Os, in case of Cr based compounds 
therefore arises due to the relative positioning of the Cr 
and B' energy levels, which narrows down the width of 
the B' states substantially in the up spin channel. The 
development of the intrinsic moment is thus helped by 
the delicate energy level structure responsible for the HD 
mechanism and would not have been present otherwise. 



moment, due to the uncompensated orbital moment at 
B' site. Whether SOC has any infiuence on the trend 
within the T^'s, therefore needs to be explored. Interplay 
of two driving mechanisms in W-Re-Os series, however, 
makes it difficult to extract the magnetic exchanges as 
energy difference between two specific magnetic configu- 
rations. The presence of finite, intrinsic moment at B' 
site, as is the case for Re and Os compounds, makes 
the moment at B' site frustrated in an antiferromagnetic 
(AFM) configuration of Cr spins. The finite presence 
of HD mechanism, on the other hand, disfavors stabi- 
lization of magnetic configurations with majority of B 
spins aligned parallely to that of the moment at B' site. 
Admitting these difficulties, we carried out total energy 
calculations of the two possible spin configurations, one 
FM arrangement and another A-type AFM arrangement 
of Cr spins, with Cr spins between two adjacent planes 
are antiferromagnetically coupled, and are ferromagnet- 
ically coupled in plane. For the AFM calculations, the 
moment at B' site was found to be antiparallely (par- 
allely) aligned to the spins of in-plane (out-of-plane) Cr 
sites which are four (two) in number. As expected, the 
energy differences between FM and AFM configurations 
are found to be positive for all cases proving FM arrange- 
ment of Cr spins to be the stable phase. The values of 
the energy difference is found to increase from W to Re 
system (from 0.23 eV/formula unit to 0.25 eV/formula 
unit), and then decrease from Re to Os system (from 
0.25 eV/formula unit to 0.24 eV/formula unit). Intro- 
duction of SO interaction though changes the individual 
energy differences by about 0.03 eV, the trend remains 
unaltered. This indicates that although the presence of 
substantial SOC at B' site is important for producing the 
large magneto-optical signals ,1^-^ it plays little role in set- 
ting up the trend in Tc- In order to examine the role of 
missing correlation effect in GGA, we have also carried 
out LDA+U calculations with a choice of U value of 3 eV 
at Cr site and 0.8 eV at B' site. Application of larger U at 
Cr site and relatively smaller one at B' sites, as expected 
from the relative band- widths, were found to preserve the 
general conclusions intact with more localized character 
of d states at B' site, as expected. 



D. Exact Diagonalization Study of Model 
Hamiltonian 



C. Total Energy Calculations 

The calculations discussed so far, do not include SOC, 
which may be important. The numbers within the 
bracket in right panel of Fig.l show the individual spin 
and orbital moments as well as net moments as obtained 
in GGA+SO calculations. As expected, the orbital mo- 
ments are large at B' site. Interestingly, we find while 
the net moment was zero for SCOO without SOC, it is 
the consideration of SOC that gives rise to a non-zero 



In view of the difficulty in stabilizing the appropri- 
ate excited state magnetic configuration, the stability of 
the FM arrangements of Cr spins may be measured as 
the energy difference between the FM spin configuration 
and the paramagnetic (PM) phase. The description of 
the PM phase needs consideration of different disordered 
spin configurations and averaging over a large number of 
them, which is almost impossible within the DFT frame- 
work. Such calculations are much easier to handle within 
a model Hamiltonian description. The model Hamilto- 
nian, describing the interplay of the HD and SE mecha- 
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nism may be written as, 



H = eCrY^ fLjtaa + fB' ^La^^'^o 



ieB 



ieB' 



i<T,a"''J<y,a 



-tec Y fL,afj<y,a + J ^ • fl^aapfip 

<ij>{T,a i^Cr 

+J2 Y ■ ^1 

ieCrJeB 

where the /'s and m's refer to the Cr t2g and B' t2g de- 
grees of freedoms. tcB', tB'B', tec represent the nearest 
neighbor Cr-B', second nearest neighbor B'-B' and Cr- 
Cr hoppings respectively, a is the spin index and a is 
the orbital index that spans the t2g manifold. The dif- 
ference between the ionic levels, A = tcr — ^B', defines 
the on-site energy difference between Cr t2g and B' t2g 
levels. Sj is the intrinsic moment at the B' site. The first, 
six terms of the Hamiltonian, represent the HD mecha- 
nism, which consist of a large core spin at the Cr site {Si) 
and the coupling between the core spin and the itinerant 
electron delocalized over the Cr-B' network. Variants of 
this part has been considered by several authorJi^Hini [n 
the context of SFMO. The last term represents the SE 
mechanism, that consists of coupling between Cr spin and 
the intrinsic moment at B' site. The parameters of the 
model Hamiltonian are extracted out of DFT calculations 
through NMTO downfolding technique of constructing 
the real space Hamiltonian in the basis of effective Cr 
t2g and B' t2g degrees of freedom, tcs', tB'B' and tec 
hoppings are found to -0.35 eV, -0.12 eV and -0.08 eV 
respectively, with little variation within the W-Re-Os se- 
ries. A-s show a varying trend within the W-Re-Os series 
(A'^ = -0.66 eV, A^'^ = 0.03 eV, A^" = 0.26 eV). The 
parameters involving J and J2 were obtained from the 
spin-splitting at Cr site and the extra splitting observed 
at B site as compared to that expected from the electron 
filling effect and the splitting for W compound. 

The constructed model is then solved using exact di- 
agonalization on a lattice of dimension 8x8x8. Cal- 
culations have been carried out as well for lattices of size 
4x4x4 and 6x6x6. The trend is found to be 
the same as presented for 8 x 8 x 8. Exact diagonal- 
ization was first carried out considering the B' site to be 
totally nonmagnetic, i.e. setting the last term to zero, 
which boils down to the same underlying model Hamilto- 
nian as that of SFMO. The energy difference between the 
PM and FM is found to decrease with increasing num- 
ber of valence electron, as shown by diamond symbols 
in Fig. 3. This is exactly the similar tre nd as found in 
a recent calculation on La doped SFM0,'2^J^ as well as 
in Ref. ,9, with FM getting destabilized with increase of 
valence electrons. This trend of suppression of upon 
increasing valence electron count is further amplified due 
to change in A within the Cr-B' series. This variation 
in A in the present series is in contrast to the prescrip- 
tion given in reflHI to achieve high Tc. Upon reaching 



Figure 3: PM - FM energy differences plotted as a function 
of valence electron count, as obtained in exact diagonaliza- 
tion calculation. Tlie diamond (square) symbols connected 
by dashed (solid) line correspond to calculations correspond- 
ing to Hamiltonian, without (with) J2 term. 



valence electron count equal to 3 which corresponds to 
Os compound, FM becomes totally unstable, refiected in 
negative sign of the energy difference. We note a rather 
rapid decrease in moving from N=2 case to N=3 case. 
This may reflect the special situation of Os compound, 
with Cr t2g-H' t2g configuration, an ideal super-exchange 
situation with insulating solution, that adds on to the 
general trend. The situation gets dramatically changed 
upon inclusion of the growing localized magnetic nature 
of the B' site, as shown by square symbols in Fig. 3. Con- 
sidering J2 values, as obtained in GGA calculations, we 
find that PM and FM energy difference, recovers the cor- 
rect trend in moving from W to Re to Os as has been ob- 
served experimentally. Mapping the PM and FM energy 
difference to the mean field Tc, one obtains values 870 
K, 1160 K and 1450 K for the W, Re and Os compounds 
respectively. Although the values are overestimated com- 
pared to experimental values, presumably due to the fi- 
nite size effect of exact diagonalization calculation and 
the mean field formula, the trend is very well reproduced 
with TfVT^ = 1.33 and T^VTf ^ = 1.25, compared to 
experimental estimates of 1.38 and 1.17 respectively.'^ts' 



IV. CONCLUSION AND DISCUSSION 

In conclusion, we have studied the counter-intuitive 
Tc trend in Cr based double perovskites, Sr2CrB'06 
(B'=W/Re/Os). Analysis of the electronic and magnetic 
properties shows that the progressive enhancement of the 
Tc across the 5d series should be understood as the inter- 
play of two driving mechanisms: HD mechanism respon- 
sible for the negative spin splitting at B' site as in SFMO, 
and SE mechanism. The HD mechanism gets weaker as 
one moves along the series from W to Re to Os, due to 
the increased energy level separation of B' from Cr. SE, 
on the other hand, gets stronger in moving from SCWO 
to SCRO to SCOO due to the presence of growing in- 
trinsic moment at B' site, following the dehybridization 
effect. The observation of uncompensated moment in 
SCOO arises due to the presence of SO. SCOO, in that 
sense, should be thought as a ferrimagnet rather than a 
ferromagnet. With this, we demystify the puzzling Tc 
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trend in Cr-B' double perovskite series. knowledges discussion with L. AlfF. 
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